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Abstract 

We explored vibrations of a single-degree of freedom oscillator with a magneto- 
rheological damper subjected to kinematic excitations. Using fast and slow scales 
decoupling procedure we derived an effective damping coefficient in the limit of 
high frequency excitation. Damping characteristics, as functions of velocity, change 
considerably especially by terminating the singular non-smoothness points. This 
effect was more transparent for a larger control parameter which was defined as the 
product of the excitation amplitude and its frequency. 
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1 Introduction 



Magneto-rheological (MR) semi-active dampers with multiple ratio charac- 
teristics are able to be controlled efficiently by applied magnetic field vibra- 
tions of various dynamical systems [1,2,3,4]. Having relatively small switching 
time (between different damping ratio) they have replaced traditional electro- 
hydraulic dampers [5]. In this paper we will focus on an effective damping 
term in the presence of fast kinematic oscillations. To describe the damp- 
ing phenomenon of the suspension quarter-car model, we use the Bingham 
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viscoplastic model [1,2], which assumes a mixture of viscous and dry friction 
terms (composed of piece linear regions). To obtain the effective damping term 
we will use the two time scales splitting which is relevant in large frequency 
kinematic excitation due to the motion of car on the rough road surface. 

The present paper in divided into four sections. Following the Introduction 
(Sec. 1) we will present the Model (Sec. 2). The calculation schema and the 
results obtained will be offered in Sec. 3, Calculations and Results. We will 
finish with Sec. 4, Summary and Conclusions. 




Fig. 1. A single degree of freedom of quarter-car model. 
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Fig. 2. Damping term K against velocity for a = 0.5 and fa = 0.1, k 2 = 0.5, 
Ak = 0.0 for '1'; fa = 0.1, k 2 = 0.5, Ak = 0.3 for '2'; fa = 0.1, k 2 = 0.0, Ak = 0.0 
for '3'; fa = 0.0, k 2 = 0.0, Ak = 0.0 for '4'. 

2 The model 

The dimensionless equation of motion for a single degree of freedom oscillator 
frequently used in description of a quarter-car model [6] 

x + K(x) + x = to 2 xq sin cut, (1) 
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where —K(x) is a damping force and u 2 xo sin cut is a kinematic excitation 
term dependent on the vehicle speed and the roughness of the road surface. 
For simplicity we assumed that this surface can be described by a harmonic 
function (Eq. 1). 
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Fig. 3. Averaged damping term < K > against velocity v where k% = 0.1, k 2 = 0.5, 
Ak = 0.0 for '1'; k x = 0.1, k 2 = 0.5, Ak = 0.3 for '2'; k x = 0.1, k 2 = 0.0, Ak = 0.0 
for '3' (as in Fig. 1), and loxq = 10. 

Here we assumed that the nonlinear damping force —K(x) is of Bingham type 
[1,2] 

K(x) = k\(l + asgnx)x + k 2 sgnx + Ak. (2) 

Some examples of possible damping characteristics are presented in Fig. 2. 
Note, the above force term possesses singular points where the damping is 
non-continuous (lines '1','2') or non-smooth (line '3'). Usually, the mirror sym- 
metry (v — > —v) is broken by the change in slopes (line '3') [6]. Additionally, it 
can be also associated with constant term ('2' and '3') related to the direction 
dependent dry friction phenomenon [7]. Using the above model we will inves- 
tigate the effect of oscillations with high frequency. For a quarter-car model 
(Fig. 1) this effect could appear during its motion with high speed. 



3 Calculations and Results 

The large value of excitation frequency uj » 1 introduces naturally to second 
'fast' time scale in our system (Fig. 1). 

In fact the multiple time-scale method is an efficient and widely used treat- 
ment for nonlinear mechanical systems [8,9,10] to get an approximate analytic 
solution. In these cases the time-scales are introduced to account for nonlinear 
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terms in the examined systems. On the other hand, the fast time scale can be 
introduced through excitation terms in the limit of very high excitation fre- 
quency [7,11,12,13]. In various nonlinear systems [7], subjected to additional 
excitation fast and slow time-scales, such splitting is realized physically. It 
enables to estimate a behaviour in terms of the averaged system equations of 
motion after averaging over fast oscillations. 

Note, comparing to our system (Eqs. 1,2) in the original treatment by Thom- 
sen, k\ = k 2 ^ and A A; ^ [7,11] to describe the effect of dry friction. On 
the other hand, in Ref. [6] ki = and Ak = 0, k\ ^ for a quarter-car model. 
Below we will consider the most general case with damping (Eq. 2) described 
by all nonzero terms. 

Assuming that u » 1 we define second (fast) time scale by introducing a 
small parameter: 

e = uj- 1 (3) 



Following Thomsen (2005) [11] we define two time scales fast given by T and 
slow by Ti, respectively 

T = e-H, T x = t. (4) 



Then two corresponding time derivatives D and D x are 

D = d/dT , D l = d/dT l . (5) 

The original time derivatives can be written as follows 

± = e~ 1 D + D 1 , (6) 
^ = e- 2 D 2 + 2e- 1 D D 1 + D 2 1 , 
while a two component solution reads 

x(T ,T 1 ) = z(T 1 ) + eZ(T ,T 1 ). (7) 



By introducing the above quantities (Eqs. 4,7) into the initial equation of 
motion (Eq. 1) and collecting terms by e _1 we obtain equation for a fast 
motion part 

DlZ = wx dn(T ). (8) 
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On the other hand for a slow motion one gets (collecting terms by e°) 

D\z + K{D x z + DoO + z = 0. (9) 



To find a solution for z in the last equation (Eq. 9) we have to examine fast 
oscillations in £ and perform averaging of the whole equation in the excitation 
period 2ir (in the fast T time scale). 

At first, we get a solution for £ from Eq. 8 




Fig. 4. Time histories of damping term for system parameters: k\ = 0.05, hi = 0, 
Ak = 0, u = 100 and (a) a = 0.1, (b) a = 0.5, (c) a = 0.1, (d) a = 0.5, and 
xq = 0.01 (red line), xo = 0.1 (green line) respectively. 

£(T ) = -^x sin(T ). (10) 
Then averaging procedure for a damping term can be written as 



< K(D lZ + DqO > = — [ + asgn(D 1 z + D ^)){D lZ + D £) 

o 

+ k 2 sgn(D lZ + £>oO) + Ak ) dr ( n ) 
By using Eq. 10 and evaluating the above integral we get 
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Fig. 5. Averaged values of a damping term < K > against velocity v (a = 0.5, 
k\ = 0.1) without excitation '1' and with an excitation term '2 '-'3' for ooxq = 5 and 
10, respectively. Note the curve '3' may correspond to case xo = 0.1 and u = 100. 

1 2n 

< K > = kia— / (Diz — x uj cos r)sgn(D 1 z — x u cos r)dr 

Z7T J 


1 

+ k 2 — / sgn(DiZ — x u;cosr)dr + Ak, (12) 
o 

where r = T G [0, 2tt]. 

After integration we get an analytical formula for the averaged damping term 



< K > 



, 2 / D 1 z 
1 arccos 

7T \X UJ 



2 

H — sin 

7T 



arccos 



/ -Pig 



(A;2 + kicxDxz) 
kiax uj + k\D\z + Afe. 



(13) 



It is obvious that Xquo is a control parameter which scaling influences averaged 
damping in a nontrivial way. 

Following Thomsen [7,11] discussion on dry friction, we decided to plot the 
results (Eq. 13) for three values of xquj and k% = and Ak = in Fig. 3. For 
comparison we have also performed the simulations of the system, using Eq. 
1, for large u (u = 100). Figures 4a-d show the results for different choice of 
system parameters. The simulation results show that the product Xqoj and a 
determine the effective damping. It is visible directly in envelops of quenching 
curves as well as in the asymptotic equilibrium value < x >= lim^ooX. In 
our simulation cases Figs. 4 a-d we have got < x >~ 0, -0.02, -0.006, - 
0.031 for red lines; < x >~ 0, -0.2, -0.06, -0.31 for green lines. Note that the 
simulation result shown in Fig. 4d < x >~ —0.31 (for the green line) agrees 
with the analytic result (Eq. 12, Fig. 3 - line '3'). This agreement is due to 
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the equilibrium of effective damping for v — and the static displacement 
force (note, the stiffness coefficient is equal to 1). Finally, in Fig. 5 we show 
the effective damping < K > for various choices of damping terms (Eq. 2) 
including asymmetric friction phenomenon related to curves plotted in Fig. 2. 
When comparing these two figures (Figs. 2 and 5) one can observe that fast 
oscillations can easily eliminate the non-smoothness and discontinuities in the 
original characteristics. As we expected, asymmetries in the initial damping 
terms (Eq. 2) lead to new equilibrium point. 



4 Summary and Conclusions 

We have examined vibrations of a single-degree of freedom oscillator with MR 
fluid damper subjected to kinematic excitations. Nontrivial behaviour of the 
system is immanent feature of the MR fluid damper. For defined magnetic 
fluid characteristics [2,6], we have used the procedure based on fast and slow 
scales splitting and, in this way, we have derived effective damping coefficient 
depending on magnetic fluid characteristics. Consequently, we have shown how 
the initial oscillations were quenched out in presence of fast kinematic forcing. 
Presence of asymmetry in the original damping terms (Eq. 2) changes the 
equilibrium point. 

We should add that the effect of fast excitations can appear in a different 
way including external forcing of various kinds and an additional degree of 
freedom. Particularly this effect can be used to improve brake systems in cars 
[14]. 
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